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Recent measurements of cosmic ray proton and helium spectra show a hardening above a few 
hundreds of GeV. This excess is hard to understand in the framework of the conventional models of 
Galactic cosmic ray production and propagation. We propose here to explain this anomaly by the 
presence of local sources. Cosmic ray propagation is described as a diffusion process taking place 
inside a two-zone magnetic halo. We calculate the proton and helium fluxes at the Earth between 
50 GeV and 100 TeV. Improving over a similar analysis, we consistently derive these fluxes by taking 
into account both local and remote sources for which a unique injection rate is assumed. We find 
cosmic ray propagation parameters for which the proton and helium spectra remarkably agree with 
the PAMELA and CREAM measurements over four decades in energy. 
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Understanding the origin and propagation of cosmic 
rays (CR) in the Galaxy is a long standing puzzle whose 
solution indeed needs the combination of several different 
observations in a wide energy range. The most realistic 
description of CR propagation is given by diffusion mod- 
els. Two main approaches have been developed so far: 
analytical (or semi-analytical) diffusion models and fully 
numerical diffusion models. Well-known realizations of 
these two different approaches are respectively the two- 
zone model with USINE code [TH3], and the GALPROP 
package HE] which has led to the DRAGON code 0[7]. 
In the conventional scheme, primary cosmic ray nuclei are 
accelerated by the shocks of the supernova explosions of 
massive stars. They are injected inside the Galactic mag- 
netic halo with a rate q oc lZ~ a where 1Z = p/Ze stands 
for the rigidity and a ~ 2.15 ± 0.15. The particles are 
subsequently scattered by the turbulent irregularities of 
the Galactic magnetic field. Their transport is described 
phenomenologically by space diffusion with a coefficient 
K oc 1Z S whose energy dependence is characterized by 
the index 8. The boron to carbon (B/C) ratio is a tracer 
of cosmic ray propagation and points towards a value of 
S £ [0.4, 0.85] for the diffusion index. At high energy, the 
flux of a given primary CR species at the Earth is given 
by $ oc q/K and falls with energy typically like 1/E a+S . 
Its spectrum exhibits a power-law behavior whose index 
is the sum of a and 5. So far, it has been generally un- 
derstood that the spectrum of primary cosmic rays ap- 
proximately behaves as E~ 2 - 7 until the knee, which is a 
bend downwards at around 3 PeV. The spectrum beyond 
the knee falls approximately like E~ 3A . 

Recent measurements of the absolute high-energy 
CR proton and helium spectra have been reported by 



CREAM E] and PAMELA [10] experiments. Obser- 
vations point towards an excess in the CR proton and 
helium fluxes above 250 GeV/nuc. The single power- 
law hypothesis is rejected at 95% C.L. The hardening 
of the proton spectrum occurs at 232^^ GeV with a 
change of the spectral index from 2.85±0.015± 0.004 to 
2.67 ± 0.03 ± 0.05. For the helium data, the spectral in- 
dex varies from 2.766± 0.01 ± 0.027 to 2.477± 0.06± 0.03 
with the hardening setting in at 2431^ GeV/nuc. These 
results challenge the conventional scenarios proposed so 
far to model Galactic cosmic rays. They indicate the 
presence of an anomalous behavior of the CR proton and 
helium spectra in the 100 GeV to 100 TeV energy range. 

Explanations of this anomaly have been tentatively 
given since its discovery. They mainly imply a mod- 
ification of the energy behavior of either the injection 
spectrum q(E) or the diffusion coefficient K(E). To 
commence, a break in a could arise from a modification 
of the conventional diffusive shock acceleration (DSA) 
scheme, as suggested by Malkov et al. |TT] and Ohira 
and Ioka [12] . In the same vein, the possibility of differ- 
ent classes of CR sources has been proposed some time 
ago by Stanev et al. [T3] and Zatsepin and Sokolskaya 
[14) . For instance, cosmic rays accelerated in the magne- 
tized winds of exploding Wolf-Rayet and red supergiant 
stars could have a double spectrum, with a hard compo- 
nent produced in the polar cap regions of these objects. 
According to Bicrmann et al. [15) . this hard component 
would take over the smooth one above a few hundreds 
of GeV, hence the observed break in the proton and he- 
lium spectra. Not so different is the proposition of Yuan 
et al. |16j where a spread in the injection index a is intro- 
duced. Another direction implies a modification of the 



2 



diffusion coefficient K. As proposed by Ave et ai. [T7] . 
the proton and helium anomaly could be due to a wel- 
come, but unexpected, decrease of the spectral index 8 
at high energy. Recently Blasi et al. [T5] have given some 
theoretical motivations to such changes in diffusion. A 
local variation of K could also have a similar effect as sug- 
gested by Tomassetti [H|. Finally, inspired by Hdrandel 
et al. 20J, Blasi and Amato [21] have invoked an unusu- 
ally strong spallation of the CR species on the Galactic 
gas. This possibility has been recently criticized in a de- 
tailed analysis carried out by Vladimirov et al. [22] of 
some of the above mentioned solutions to the CR proton 
and helium anomaly. 

In this Letter, we show that the proton and helium 
spectral hardening above 250 GeV/nuc can be attributed 
to local sources of cosmic rays, whose presence is associ- 
ated to known supernova remnants (SNR) and pulsars. 
These objects can be found in astronomical catalogs such 
as the Green catalog [23] which can be completed with 
the ATNF pulsar database [211 [25] . In our approach, 
there is no need to modify the conventional CR prop- 
agation model. In particular, the variations with CR 
energy of the injection rate q of individual sources and 
of the space diffusion coefficient K are power laws re- 
spectively characterized by the spectral indices a and S. 
This idea has already been suggested recently by Erlykin 
and Wolfendale |26| who explain the hardening with very 
few sources (mainly Monogem Ring) and by Thoudam 
and Hdrandel [27J who consider a catalog of 10 nearby 
sources. The principal weakness of these analyses is the 
lack of a consistent treatment of the CR spectra in the 
entire energy range extending from tens of GeV up to 
a few PeV. This is particularly clear in [27] where the 
proton and helium anomaly is derived from a handful of 
local sources, whereas the low energy spectra of these 
species are not calculated but merely fitted in order to 
get a value for a once 5 has been chosen. It should be 
noted that the magnitude of the CR proton (helium) flux 
is related over the entire energy range to the injection 
rate q of individual sources. The low energy (power-law 
regime) and high energy (spectral hardening) parts of the 
CR spectra are connected with each other. A consistent 
treatment of the problem requires that the proton and 
helium fluxes are calculated over the entire energy range. 
A crucial problem is also to understand why just a few 
local sources could explain the spectral hardening at high 
energies whereas the bulk of the Galactic sources is re- 
quired in order to account for the power-law behavior of 
the fluxes below 250 GeV/nuc. This aspect, which is not 
addressed in the above mentioned analyses, bears upon 
the more general question of the discreet nature of the 
sources. In the conventional model of CR propagation, 
these are treated as a jelly spreading over the Galac- 
tic disk and continuously accelerating cosmic rays. The 
question arises then to understand why and in which con- 
ditions that scheme breaks down at high energies where 



local and point-like objects come into play. The results 
presented here are based on a detailed investigation of 
that question. Bernard et al. [28j have recently shown 
how to reconcile the presence of point-like sources with 
the conventional description of CR production and prop- 
agation. We briefly recall the salient features of their 
analysis. 

Once accelerated by the sources that lie within the 
Galactic disk, CR nuclei diffuse on the irregularities of 
the Galactic magnetic field. The diffusion coefficient 
K = K fllZ s accounts for that process, where K is a 
normalization constant and fl denotes the particle veloc- 
ity. The magnetic halo, inside which cosmic rays propa- 
gate before escaping into intergalactic space, is assumed 
to be a flat cylindrical domain which matches the circular 
structure of the Milky Way. The Galactic disk is sand- 
wiched between two confinement layers whose thickness 
L is unknown and turns out to be crucial in our inves- 
tigation. Stellar winds combine to generate a Galactic 
convection that wipes cosmic rays away from the disk, 
with velocity V c (z) = V c sign(,z). CR nuclei also un- 
dergo collisions with the interstellar medium (ISM) with 
a rate T sp = a co \ fl tiism ■ Above a few GeV, diffusive re- 
acceleration and energy losses may be disregarded and 
the master equation for the space and energy number 
density ij) = dn/dT of a given CR species simplifies into 
the diffusion equation 

^ + d z (V c t/,) - K(E) + r Bp V = q acc . (1) 

The CR transport parameters K 0l 6, L and V c can be 
weakly constrained from the B/C ratio (see for instance 
[TJ [2S] [3D] ) • In the conventional approach, the CR source 
term q acc is a continuous function of space and time. 
Steady-state is moreover assumed. This is an oversim- 
plification insofar as CR sources are actually point-like, 
with an average supernova explosion rate v of 1 to 3 
events per century. In the stochastic treatment devel- 
oped by Bernard et al. [25], the production rate of CR 
nuclei through acceleration is given by 

9acc(x S ,ts) = ^2 Qi S 3 (*S ~*i)$(ts -U) , (2) 

lev 

where each source i that belongs to the population V 
contributes a factor qi at position Xj and time U. The 
total flux <3? = (l/47r) flip at the Earth depends on the 
precise locations and ages of all the sources and varies 
from one particular population V to another. Because 
we do not know the actual distribution of the Galactic 
sources that have generated the observed CR flux, we 
must rely on a statistical analysis and consider the po- 
sition and age of each source as random variables. The 
CR flux &(E) at a given energy E behaves as a stochas- 
tic variable whose probability distribution function p($) 
has been studied in [25]. The conventional CR model 
is recovered by taking the statistical average of the flux 
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FIG. 1: Proton (left panel) and helium (right panel) spectra in the range extending from 5 GeV/nuc to 100 TeV/nuc. The 
solid purple line 4> is the mean flux produced by a continuous distribution of all the Galactic sources. It corresponds to the 
prediction of the conventional steady-state propagation model. The dotted short-dashed green curve $i oc is the mean flux 
obtained by randomly picking local sources (less than 2 kpc away and younger than 30,000 years). The long-dashed red line 
"I'cat is produced by the sources of our neighborhood found in SNR and pulsars catalogs. Twenty objects contribute to the 
signal displayed here. The dotted long-dashed black curve features the sum of the contributions from the local catalog sources 
<E>cat and from the continuous distribution of the far away sources $ C xt- Solar modulation is switched on with a Fisk potential 
4>y of 475 MV. Data are from PAMELA [10] and CREAM [8]. The CR propagation parameters used to get this plot correspond 
to model A of Tab. |U 



over the ensemble of all possible populations V . This av- 
erage flux $ turns out to be the solution of Eq. [l] with a 
continuous source term t7 acc . More exciting is the spread 
of the flux $ around its average value $. Using a Monte 
Carlo approach, Bernard et al. [28] have shown that if the 
magnetic halo is thin, the statistical fluctuations of the 
flux may be significant. The residence time Tdn ~ L 2 / K 
of CR nuclei within the magnetic halo decreases with its 
thickness L. The number N of sources that contribute 
to the signal at the Earth scales as vr^f. The smaller L, 
the smaller N and the larger the flux variance. Should 
the magnetic halo be sufficiently thin, we expect fluctu- 
ations of the flux, especially at high energies where K 
becomes large. In this case, the hardening of the proton 
and helium spectra appears to be a mere fluctuation of 
the CR flux whose probability to occur is not vanishingly 
small. According to this line of reasoning, the proton and 
helium anomaly results from the particular configuration 
of the actual CR sources. These objects are incidentally 
known in the nearby region for which catalogs of SNR 
and pulsars are available. The domain extending 2 kpc 
around the Earth and encompassing objects that have ex- 
ploded less than 30,000 years ago is defined as the local 
region. The catalogs are no longer complete outside and 
fail to be reliable. In the conventional CR model, the 



local sources would yield an average contribution <&i oc 
whereas the actual objects yield a much larger flux <& cat . 
Denoting by $ cxt the flux from the other sources, we infer 
a total signal at the Earth 

* = $cat + *cxt , (3) 

to be compared to the prediction of the conventional 
steady-state model 

i = i>i oc + $ oxt . (4) 

The flux produced by the external sources has a very 
small variance as shown in |28j . We may then identify 
$ cxt with its statistical average <& ex t- 

We have performed a scan over the CR parameters 
in order to fit the PAMELA [10] and CREAM [8] data 
over an energy range extending from 50 GeV/nuc to 
100 TeV/nuc. For each model, CR propagation is de- 
scribed by Kq, S and L. At high energy, Galactic con- 
vection and solar modulation have little effect on the CR 
flux. That is why we have set V c and the Fisk poten- 
tial </>f equal to during the scan. The injection indices 
a p and a^ c are independently adjusted to the proton 
and helium spectra. Observations point towards slightly 
different power laws for these fluxes at energies below 
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CR model 


K [kp^Myr- 1 ] 
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L [kpc] 


V c [km s" 1 ] 


a P 




v [century ] 


xl 




A 


1.29 x 10" 2 


0.63 


1 





2.23 


2.14 


1.09 


0.33 


0.54 


B 


5.10 x 10" 3 


0.75 


2 


11.5 


2.11 


2.0 


0.8 


0.28 


0.59 



TABLE I: The two sets of CR injection and propagation parameters featured in this table provide very good fits to the 
PAMELA and CREAM data from 50 GeV/nuc to 100 TeV/nuc. The proton and helium fluxes are simultaneously adjusted 
with the same values of Ko, S and L. The injection indices a p and «h c are determined independently from each other. The 
average supernova explosion rate is denoted by v. The results of the fits to the proton and helium spectra are respectively 
gauged by x% an d Xh c - Model A corresponds to the best fit and is featured in Fig. [I] The agreement with the observations is 
remarkable even down to 1 GeV and spans five orders of magnitude. Model B belongs to a set of CR parameters which have 
been shown in [I] to be compatible with the B/C ratio. Although the total \ 2 is the same as for model A, the agreement with 
the data is not as good below 100 GeV. 



250 GeV/nuc. The last injection parameter which we 
have considered is the average supernova explosion rate 
v in the Galaxy. Each CR configuration is then char- 
acterized by six parameters. The quality of the fits to 
the proton and helium data is respectively gauged by the 
reduced chi-squares \\ an d Xhc- We have found many 
models which reproduce fairly well the power-law regime 
and the hardening of both the proton and helium fluxes. 
The two configurations of Tab. [I] adjust particularly well 
the PAMELA and CREAM spectra. 

This excellent agreement makes us confident that the 
proton and helium anomaly can actually be explained 
by existing local sources which have been extracted from 
SNR and pulsar surveys. The model which we have pre- 
sented here is quite simple. Refining it is beyond the 
scope of this Letter. Some directions can nevertheless 
be given in order to improve the solution which we have 
just sketched. To commence, the best fit is obtained for 
a rather small value of the magnetic halo thickness L. 
This is not only true for models A and B but is a gen- 
eral trend easy to understand. As already explained, the 
thinner the magnetic halo, the smaller the number N of 
sources which contribute to the total signal and the larger 
the injection rate q of individual sources. The contribu- 
tions $ cat and $i oc from the local region are no longer 
swamped in the total flux when L is small. This may 
be a problem as recent studies [22 [22] of the gamma-ray 
and synchrotron diffuse emissions seem to favour rather 
large values of L. A possible improvement of our model 
would be to distribute the CR sources within spiral arms 
and to take into account the rotation of the Galaxy. The 
sources which we have considered here in order to derive 
the contribution $ oxt to the flux are equally spread along 
the azimuthal direction. It would be interesting to inves- 
tigate if $oxt decreases in a more realistic setup. Notice 
also that in order to get a significant injection rate q, we 
are naturally driven towards a small supernova explosion 
rate v. The values found for models A and B are close 
to 1 explosion per century, at the lower edge of the plau- 
sible range. The local sources which we have extracted 
from the catalogs correspond to a larger rate v of 3.3 
events per century. Because the Sun lies near two Galac- 
tic arms, the average explosion rate in our neighborhood 



could reasonably be higher than the mean rate of the 
Galaxy. This is actually supported by our catalog. As 
shown in Fig. 7 of |28j , the number of known sources in 
our vicinity is compatible with a value of v larger than 3 
explosions per century for the past 3 x 10 4 years (depend- 
ing on the radial distribution of the CR sources along the 
Galactic disk). Taking into account the Galactic spiral 
arms and their rotation could lead to a larger value of the 
average explosion rate and alleviate the apparent discrep- 
ancy between the average and local values of v. Finally, 
we have modeled the supernova explosions as point-like 
events. Cosmic rays are believed to be accelerated in the 
shocks which follow these explosions and which propa- 
gate in the ISM during 10 5 years. The injection sites are 
more spherical shells than points. Depending on the CR 
energy, Thoudam and Horandel [57] quote escape times 
between 500 and 10 5 years after the stellar explosion. 
The injection takes place from a remnant whose radius 
varies from 5 to 100 pc. Taking into account the actual 
structure of CR accelerators could substantially modify 
the contribution $ cat to the total signal, allowing larger 
values of L to provide acceptable fits to the PAMELA 
and CREAM data. The simplistic solution to the proton 
and helium anomaly which we have sketched in this Let- 
ter is definitely exciting in spite of the above mentioned 
problems and should motivate further investigations. 
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